Clinical associations with a polygenic predisposition to benign lower white blood cell counts.

Polygenic variation unrelated to disease contributes to interindividual variation in baseline white blood cell (WBC) counts, but its clinical significance is uncharacterized. We investigated the clinical consequences of a genetic predisposition toward lower WBC counts among 89,559 biobank participants from tertiary care centers using a polygenic score for WBC count (PGSWBC) comprising single nucleotide polymorphisms not associated with disease. A predisposition to lower WBC counts was associated with a decreased risk of identifying pathology on a bone marrow biopsy performed for a low WBC count (odds-ratio = 0.55 per standard deviation increase in PGSWBC [95%CI, 0.30-0.94], p = 0.04), an increased risk of leukopenia (a low WBC count) when treated with a chemotherapeutic (n = 1724, hazard ratio [HR] = 0.78 [0.69-0.88], p = 4.0 × 10-5) or immunosuppressant (n = 354, HR = 0.61 [0.38-0.99], p = 0.04). A predisposition to benign lower WBC counts was associated with an increased risk of discontinuing azathioprine treatment (n = 1,466, HR = 0.62 [0.44-0.87], p = 0.006). Collectively, these findings suggest that there are genetically predisposed individuals who are susceptible to escalations or alterations in clinical care that may be harmful or of little benefit.

Clinical consequences of a polygenic predisposition to benign lower white blood cell counts: Consequences of benign WBC count genetics.

Polygenic variation unrelated to disease contributes to interindividual variation in baseline white blood cell (WBC) counts, but its clinical significance is undefined. We investigated the clinical consequences of a genetic predisposition toward lower WBC counts among 89,559 biobank participants from tertiary care centers using a polygenic score for WBC count (PGSWBC) comprising single nucleotide polymorphisms not associated with disease. A predisposition to lower WBC counts was associated with a decreased risk of identifying pathology on a bone marrow biopsy performed for a low WBC count (odds-ratio=0.55 per standard deviation increase in PGSWBC [95%CI, 0.30 - 0.94], p=0.04), an increased risk of leukopenia (a low WBC count) when treated with a chemotherapeutic (n=1,724, hazard ratio [HR]=0.78 [0.69 - 0.88], p=4.0×10-5) or immunosuppressant (n=354, HR=0.61 [0.38 - 0.99], p=0.04). A predisposition to benign lower WBC counts was associated with an increased risk of discontinuing azathioprine treatment (n=1,466, HR=0.62 [0.44 - 0.87], p=0.006). Collectively, these findings suggest that a WBC count polygenic score identifies individuals who are susceptible to escalations or alterations in clinical care that may be harmful or of little benefit.

Decomposing interaction and mediating effects of race/ethnicity and circulating blood levels of cystatin C on cognitive status in the United States health and retirement study.

Background and objectives: Elevated circulating cystatin C is associated with cognitive impairment in non-Hispanic Whites, but its role in racial disparities in dementia is understudied. In a nationally representative sample of older non-Hispanic White, non-Hispanic Black, and Hispanic adults in the United States, we use mediation-interaction analysis to understand how racial disparities in the cystatin C physiological pathway may contribute to racial disparities in prevalent dementia. Methods: In a pooled cross-sectional sample of the Health and Retirement Study (n = 9,923), we employed Poisson regression to estimate prevalence ratios and to test the relationship between elevated cystatin C (>1.24 vs. ≤1.24 mg/L) and impaired cognition, adjusted for demographics, behavioral risk factors, other biomarkers, and chronic conditions. Self-reported racialized social categories were a proxy measure for exposure to racism. We calculated additive interaction measures and conducted four-way mediation-interaction decomposition analysis to test the moderating effect of race/ethnicity and mediating effect of cystatin C on the racial disparity. Results: Overall, elevated cystatin C was associated with dementia (prevalence ratio [PR] = 1.2; 95% CI: 1.0, 1.5). Among non-Hispanic Black relative to non-Hispanic White participants, the relative excess risk due to interaction was 0.7 (95% CI: -0.1, 2.4), the attributable proportion was 0.1 (95% CI: -0.2, 0.4), and the synergy index was 1.1 (95% CI: 0.8, 1.8) in a fully adjusted model. Elevated cystatin C was estimated to account for 2% (95% CI: -0, 4%) for the racial disparity in prevalent dementia, and the interaction accounted for 8% (95% CI: -5, 22%). Analyses for Hispanic relative to non-white participants suggested moderation by race/ethnicity, but not mediation. Discussion: Elevated cystatin C was associated with dementia prevalence. Our mediation-interaction decomposition analysis suggested that the effect of elevated cystatin C on the racial disparity might be moderated by race/ethnicity, which indicates that the racialization process affects not only the distribution of circulating cystatin C across minoritized racial groups, but also the strength of association between the biomarker and dementia prevalence. These results provide evidence that cystatin C is associated with adverse brain health and this effect is larger than expected for individuals racialized as minorities had they been racialized and treated as non-Hispanic White.

Fibromuscular Dysplasia and Abdominal Aortic Aneurysms Are Dimorphic Sex-Specific Diseases With Shared Complex Genetic Architecture.

BACKGROUND: The risk of arterial diseases may be elevated among family members of individuals having multifocal fibromuscular dysplasia (FMD). We sought to investigate the risk of arterial diseases in families of individuals with FMD. METHODS: Family histories for 73 probands with FMD were obtained, which included an analysis of 463 total first-degree relatives focusing on FMD and related arterial disorders. A polygenic risk score for FMD (PRSFMD) was constructed from prior genome-wide association findings of 584 FMD cases and 7139 controls and evaluated for association with an abdominal aortic aneurysm (AAA) in a cohort of 9693 AAA cases and 294 049 controls. A previously published PRSAAA was also assessed among the FMD cases and controls. RESULTS: Of all first degree relatives of probands, 9.3% were diagnosed with FMD, aneurysms, and dissections. Aneurysmal disease occurred in 60.5% of affected relatives and 5.6% of all relatives. Among 227 female first-degree relatives of probands, 4.8% (11) had FMD, representing a relative risk (RR)FMD of 1.5 ([95% CI, 0.75-2.8]; P=0.19) compared with the estimated population prevalence of 3.3%, though not of statistical significance. Of all fathers of FMD probands, 11% had AAAs resulting in a RRAAA of 2.3 ([95% CI, 1.12-4.6]; P=0.014) compared with population estimates. The PRSFMD was found to be associated with an AAA (odds ratio, 1.03 [95% CI, 1.01-1.05]; P=2.6×10-3), and the PRSAAA was found to be associated with FMD (odds ratio, 1.53 [95% CI, 1.2-1.9]; P=9.0×10-5) as well. CONCLUSIONS: FMD and AAAs seem to be sex-dimorphic manifestations of a heritable arterial disease with a partially shared complex genetic architecture. Excess risk of having an AAA according to a family history of FMD may justify screening in family members of individuals having FMD.

Limited clinical utility for GWAS or polygenic risk score for postoperative acute kidney injury in non-cardiac surgery in European-ancestry patients.

BACKGROUND: Prior studies support a genetic basis for postoperative acute kidney injury (AKI). We conducted a genome-wide association study (GWAS), assessed the clinical utility of a polygenic risk score (PRS), and estimated the heritable component of AKI in patients who underwent noncardiac surgery. METHODS: We performed a retrospective large-scale genome-wide association study followed by a meta-analysis of patients who underwent noncardiac surgery at the Vanderbilt University Medical Center ("Vanderbilt" cohort) or Michigan Medicine, the academic medical center of the University of Michigan ("Michigan" cohort). In the Vanderbilt cohort, the relationship between polygenic risk score for estimated glomerular filtration rate and postoperative AKI was also tested to explore the predictive power of aggregating multiple common genetic variants associated with AKI risk. Similarly, in the Vanderbilt cohort genome-wide complex trait analysis was used to estimate the heritable component of AKI due to common genetic variants. RESULTS: The study population included 8248 adults in the Vanderbilt cohort (mean [SD] 58.05 [15.23] years, 50.2% men) and 5998 adults in Michigan cohort (56.24 [14.76] years, 49% men). Incident postoperative AKI events occurred in 959 patients (11.6%) and in 277 patients (4.6%), respectively. No loci met genome-wide significance in the GWAS and meta-analysis. PRS for estimated glomerular filtration rate explained a very small percentage of variance in rates of postoperative AKI and was not significantly associated with AKI (odds ratio 1.050 per 1 SD increase in polygenic risk score [95% CI, 0.971-1.134]). The estimated heritability among common variants for AKI was 4.5% (SE = 4.5%) suggesting low heritability. CONCLUSION: The findings of this study indicate that common genetic variation minimally contributes to postoperative AKI after noncardiac surgery, and likely has little clinical utility for identifying high-risk patients.

Lack of association of CYP2B6 pharmacogenetics with cyclophosphamide toxicity in patients with cancer.

PURPOSE: Cyclophosphamide is a commonly used cancer agent that is metabolically activated by polymorphic enzymes. This study aims to investigate the association between predicted activity of candidate pharmacogenes with severe toxicity during cyclophosphamide treatment. METHODS: Genome-wide genetic data was collected from an institutional genetic data repository for CYP2B6, CYP3A4, CYP2C9, CYP2C19, GSTA1, GSTP1, ALDH1A1, ALDH3A1, ABCC1, ABCB1, and ERCC1. Treatment and toxicity data were retrospectively collected from the patient's medical record. The a priori selected primary hypothesis was that patients who have CYP2B6 reduced metabolizer activity (poor or intermediate (PM/IM) vs. normal (NM) metabolizer) have lower risk of severe toxicity or cyclophosphamide treatment modification due to toxicity. RESULTS: In the primary analysis of 510 cyclophosphamide-treated patients with available genetic data, there was no difference in the odds of severe toxicity or treatment modification due to toxicity in CYP2B6 PM/IM vs. NM (odds ratio = 0.97, 95% Confidence Interval: 0.62-1.50, p = 0.88). In an exploratory, statistically uncorrected secondary analysis, carriers of the ALDH1A1 rs8187996 variant had a lower risk of the primary toxicity endpoint compared with wild-type homozygous patients (odds ratio = 0.31, 95% Confidence Interval: 0.09-0.78, p = 0.028). None of the other tested phenotypes or genotypes was associated with the primary or secondary endpoints in unadjusted analysis (all p > 0.05). CONCLUSION: The finding that patients who carry ALDH1A1 rs8187996 may have a lower risk of cyclophosphamide toxicity than wild-type patients contradicts a prior finding for this variant and should be viewed with skepticism. We found weak evidence that any of these candidate pharmacogenetic predictors of cyclophosphamide toxicity may be useful to personalize cyclophosphamide dosing to optimize therapeutic outcomes in patients with cancer.

Genetic investigation of fibromuscular dysplasia identifies risk loci and shared genetics with common cardiovascular diseases.

Fibromuscular dysplasia (FMD) is an arteriopathy associated with hypertension, stroke and myocardial infarction, affecting mostly women. We report results from the first genome-wide association meta-analysis of six studies including 1556 FMD cases and 7100 controls. We find an estimate of SNP-based heritability compatible with FMD having a polygenic basis, and report four robustly associated loci (PHACTR1, LRP1, ATP2B1, and LIMA1). Transcriptome-wide association analysis in arteries identifies one additional locus (SLC24A3). We characterize open chromatin in arterial primary cells and find that FMD associated variants are located in arterial-specific regulatory elements. Target genes are broadly involved in mechanisms related to actin cytoskeleton and intracellular calcium homeostasis, central to vascular contraction. We find significant genetic overlap between FMD and more common cardiovascular diseases and traits including blood pressure, migraine, intracranial aneurysm, and coronary artery disease.

Genome-Wide Association Study of Pelvic Organ Prolapse Using the Michigan Genomics Initiative.

OBJECTIVES: The aim of this study was to (1) replicate previously identified genetic variants significantly associated with pelvic organ prolapse and (2) identify new genetic variants associated with pelvic organ prolapse using a genome-wide association study. METHODS: Using our institution's database linking genetic and clinical data, we identified 1,329 women of European ancestry with an International Classification of Diseases, Ninth Revision (ICD-9)/ICD-10 code for prolapse, 767 of whom also had Current Procedural Terminology (CPT)/ICD-9/ICD-10 procedure codes for prolapse surgery, and 16,383 women of European ancestry older than 40 years without a prolapse diagnosis code as controls. Patients were genotyped using the Illumina HumanCoreExome chip and imputed to the Haplotype Reference Consortium. We tested 20 million single nucleotide polymorphisms (SNPs) for association with pelvic organ prolapse adjusting for relatedness, age, chip version, and 4 principal components. We compared our results with 18 previously identified genome-wide significant SNPs from the UK Biobank, Commun Biol (2020;3:129), and Obstet Gynecol (2011;118:1345-1353). RESULTS: No variants achieved genome-wide significance (P = 5 × 10-8). However, we replicated 4 SNPs with biologic plausibility at nominal significance (P ≤ 0.05): rs12325192 (P = 0.002), rs9306894 (P = 0.05), rs1920568 (P = 0.034), and rs1247943 (P = 0.041), which were all intergenic and nearest the genes SALL1, GDF7, TBX5, and TBX5, respectively. CONCLUSIONS: Our replication of 4 biologically plausible previously reported SNPs provides further evidence for a genetic contribution to prolapse, specifically that rs12325192, rs9306894, rs1920568, and rs1247943 may contribute to susceptibility for prolapse. These and previously reported associations that have not yet been replicated should be further explored in larger, more diverse cohorts, perhaps through meta-analysis.

A Novel Recurrent COL5A1 Genetic Variant Is Associated With a Dysplasia-Associated Arterial Disease Exhibiting Dissections and Fibromuscular Dysplasia.

OBJECTIVE: While rare variants in the COL5A1 gene have been associated with classical Ehlers-Danlos syndrome and rarely with arterial dissections, recurrent variants in COL5A1 underlying a systemic arteriopathy have not been described. Monogenic forms of multifocal fibromuscular dysplasia (mFMD) have not been previously defined. Approach and Results: We studied 4 independent probands with the COL5A1 pathogenic variant c.1540G>A, p.(Gly514Ser) who presented with arterial aneurysms, dissections, tortuosity, and mFMD affecting multiple arteries. Arterial medial fibroplasia and smooth muscle cell disorganization were confirmed histologically. The COL5A1 c.1540G>A variant is predicted to be pathogenic in silico and absent in gnomAD. The c.1540G>A variant is on a shared 160.1 kb haplotype with 0.4% frequency in Europeans. Furthermore, exome sequencing data from a cohort of 264 individuals with mFMD were examined for COL5A1 variants. In this mFMD cohort, COL5A1 c.1540G>A and 6 additional relatively rare COL5A1 variants predicted to be deleterious in silico were identified and were associated with arterial dissections (P=0.005). CONCLUSIONS: COL5A1 c.1540G>A is the first recurring variant recognized to be associated with arterial dissections and mFMD. This variant presents with a phenotype reminiscent of vascular Ehlers-Danlos syndrome. A shared haplotype among probands supports the existence of a common founder. Relatively rare COL5A1 genetic variants predicted to be deleterious by in silico analysis were identified in ≈2.7% of mFMD cases, and as they were enriched in patients with arterial dissections, may act as disease modifiers. Molecular testing for COL5A1 should be considered in patients with a phenotype overlapping with vascular Ehlers-Danlos syndrome and mFMD.

Chromosome 1q21.2 and additional loci influence risk of spontaneous coronary artery dissection and myocardial infarction.

Spontaneous coronary artery dissection (SCAD) is a non-atherosclerotic cause of myocardial infarction (MI), typically in young women. We undertook a genome-wide association study of SCAD (Ncases = 270/Ncontrols = 5,263) and identified and replicated an association of rs12740679 at chromosome 1q21.2 (Pdiscovery+replication = 2.19 × 10-12, OR = 1.8) influencing ADAMTSL4 expression. Meta-analysis of discovery and replication samples identified associations with P < 5 × 10-8 at chromosome 6p24.1 in PHACTR1, chromosome 12q13.3 in LRP1, and in females-only, at chromosome 21q22.11 near LINC00310. A polygenic risk score for SCAD was associated with (1) higher risk of SCAD in individuals with fibromuscular dysplasia (P = 0.021, OR = 1.82 [95% CI: 1.09-3.02]) and (2) lower risk of atherosclerotic coronary artery disease and MI in the UK Biobank (P = 1.28 × 10-17, HR = 0.91 [95% CI :0.89-0.93], for MI) and Million Veteran Program (P = 9.33 × 10-36, OR = 0.95 [95% CI: 0.94-0.96], for CAD; P = 3.35 × 10-6, OR = 0.96 [95% CI: 0.95-0.98] for MI). Here we report that SCAD-related MI and atherosclerotic MI exist at opposite ends of a genetic risk spectrum, inciting MI with disparate underlying vascular biology.

GWAS of thyroid stimulating hormone highlights pleiotropic effects and inverse association with thyroid cancer.

Thyroid stimulating hormone (TSH) is critical for normal development and metabolism. To better understand the genetic contribution to TSH levels, we conduct a GWAS meta-analysis at 22.4 million genetic markers in up to 119,715 individuals and identify 74 genome-wide significant loci for TSH, of which 28 are previously unreported. Functional experiments show that the thyroglobulin protein-altering variants P118L and G67S impact thyroglobulin secretion. Phenome-wide association analysis in the UK Biobank demonstrates the pleiotropic effects of TSH-associated variants and a polygenic score for higher TSH levels is associated with a reduced risk of thyroid cancer in the UK Biobank and three other independent studies. Two-sample Mendelian randomization using TSH index variants as instrumental variables suggests a protective effect of higher TSH levels (indicating lower thyroid function) on risk of thyroid cancer and goiter. Our findings highlight the pleiotropic effects of TSH-associated variants on thyroid function and growth of malignant and benign thyroid tumors.

Genome-wide association meta-analyses combining multiple risk phenotypes provide insights into the genetic architecture of cutaneous melanoma susceptibility.

Most genetic susceptibility to cutaneous melanoma remains to be discovered. Meta-analysis genome-wide association study (GWAS) of 36,760 cases of melanoma (67% newly genotyped) and 375,188 controls identified 54 significant (P < 5 × 10-8) loci with 68 independent single nucleotide polymorphisms. Analysis of risk estimates across geographical regions and host factors suggests the acral melanoma subtype is uniquely unrelated to pigmentation. Combining this meta-analysis with GWAS of nevus count and hair color, and transcriptome association approaches, uncovered 31 potential secondary loci for a total of 85 cutaneous melanoma susceptibility loci. These findings provide insights into cutaneous melanoma genetic architecture, reinforcing the importance of nevogenesis, pigmentation and telomere maintenance, together with identifying potential new pathways for cutaneous melanoma pathogenesis.

Exploring various polygenic risk scores for skin cancer in the phenomes of the Michigan genomics initiative and the UK Biobank with a visual catalog: PRSWeb.

Polygenic risk scores (PRS) are designed to serve as single summary measures that are easy to construct, condensing information from a large number of genetic variants associated with a disease. They have been used for stratification and prediction of disease risk. The primary focus of this paper is to demonstrate how we can combine PRS and electronic health records data to better understand the shared and unique genetic architecture and etiology of disease subtypes that may be both related and heterogeneous. PRS construction strategies often depend on the purpose of the study, the available data/summary estimates, and the underlying genetic architecture of a disease. We consider several choices for constructing a PRS using data obtained from various publicly-available sources including the UK Biobank and evaluate their abilities to predict not just the primary phenotype but also secondary phenotypes derived from electronic health records (EHR). This study was conducted using data from 30,702 unrelated, genotyped patients of recent European descent from the Michigan Genomics Initiative (MGI), a longitudinal biorepository effort within Michigan Medicine. We examine the three most common skin cancer subtypes in the USA: basal cell carcinoma, cutaneous squamous cell carcinoma, and melanoma. Using these PRS for various skin cancer subtypes, we conduct a phenome-wide association study (PheWAS) within the MGI data to evaluate PRS associations with secondary traits. PheWAS results are then replicated using population-based UK Biobank data and compared across various PRS construction methods. We develop an accompanying visual catalog called PRSweb that provides detailed PheWAS results and allows users to directly compare different PRS construction methods.

Association of Polygenic Risk Scores for Multiple Cancers in a Phenome-wide Study: Results from The Michigan Genomics Initiative.

Health systems are stewards of patient electronic health record (EHR) data with extraordinarily rich depth and breadth, reflecting thousands of diagnoses and exposures. Measures of genomic variation integrated with EHRs offer a potential strategy to accurately stratify patients for risk profiling and discover new relationships between diagnoses and genomes. The objective of this study was to evaluate whether polygenic risk scores (PRS) for common cancers are associated with multiple phenotypes in a phenome-wide association study (PheWAS) conducted in 28,260 unrelated, genotyped patients of recent European ancestry who consented to participate in the Michigan Genomics Initiative, a longitudinal biorepository effort within Michigan Medicine. PRS for 12 cancer traits were calculated using summary statistics from the NHGRI-EBI catalog. A total of 1,711 synthetic case-control studies was used for PheWAS analyses. There were 13,490 (47.7%) patients with at least one cancer diagnosis in this study sample. PRS exhibited strong association for several cancer traits they were designed for, including female breast cancer, prostate cancer, melanoma, basal cell carcinoma, squamous cell carcinoma, and thyroid cancer. Phenome-wide significant associations were observed between PRS and many non-cancer diagnoses. To differentiate PRS associations driven by the primary trait from associations arising through shared genetic risk profiles, the idea of "exclusion PRS PheWAS" was introduced. Further analysis of temporal order of the diagnoses improved our understanding of these secondary associations. This comprehensive PheWAS used PRS instead of a single variant.

Exome-wide association study reveals novel psoriasis susceptibility locus at TNFSF15 and rare protective alleles in genes contributing to type I IFN signalling.

Psoriasis is a common inflammatory skin disorder for which multiple genetic susceptibility loci have been identified, but few resolved to specific functional variants. In this study, we sought to identify common and rare psoriasis-associated gene-centric variation. Using exome arrays we genotyped four independent cohorts, totalling 11 861 psoriasis cases and 28 610 controls, aggregating the dataset through statistical meta-analysis. Single variant analysis detected a previously unreported risk locus at TNFSF15 (rs6478108; P = 1.50 × 10-8, OR = 1.10), and association of common protein-altering variants at 11 loci previously implicated in psoriasis susceptibility. We validate previous reports of protective low-frequency protein-altering variants within IFIH1 (encoding an innate antiviral receptor) and TYK2 (encoding a Janus kinase), in each case establishing a further series of protective rare variants (minor allele frequency < 0.01) via gene-wide aggregation testing (IFIH1: pburden = 2.53 × 10-7, OR = 0.707; TYK2: pburden = 6.17 × 10-4, OR = 0.744). Both genes play significant roles in type I interferon (IFN) production and signalling. Several of the protective rare and low-frequency variants in IFIH1 and TYK2 disrupt conserved protein domains, highlighting potential mechanisms through which their effect may be exerted.

A rare coding allele in IFIH1 is protective for psoriatic arthritis.

OBJECTIVES: Psoriatic arthritis (PsA) is an inflammatory arthritis associated with psoriasis. While many common risk alleles have been reported for association with PsA as well as psoriasis, few rare coding alleles have yet been identified. METHODS: To identify rare coding variation associated with PsA risk or protection, we genotyped 41 267 variants with the exome chip and investigated association within an initial cohort of 1980 PsA cases and 5913 controls. Genotype data for an independent cohort of 2234 PsA cases and 5708 controls was also made available, allowing for a meta-analysis to be performed with the discovery dataset. RESULTS: We identified an association with the rare variant rs35667974 (p=2.39x10-6, OR=0.47), encoding an Ile923Val amino acid change in the IFIH1 gene protein product. The association was reproduced in our independent cohort, which reached a high level of significance on meta-analysis with the discovery and replication datasets (p=4.67x10-10). We identified a strong association with IFIH1 when performing multiple-variant analysis (p=6.77x10-6), and found evidence of independent effects between the rare allele and the common PsA variant at the same locus. CONCLUSION: For the first time, we report a rare coding allele in IFIH1 to be protective for PsA. This rare allele has also been identified to have the same direction of effect on type I diabetes and psoriasis. While this association further supports existing evidence for IFIH1 as a causal gene for PsA, mechanistic studies will need to be pursued to confirm that IFIH1 is indeed causal.